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Abstract 

The representation of polynomials by arithmetic circuits evaluating them 
is an alternative data structure which allowed considerable progress in poly- 
CN nomial equation solving in the last fifteen years. We present a circuit based 

computation model which captures all known symbolic elimination algorithms 
^N) in effective algebraic geometry and show the intrinsically exponential complex- 

$_i ity character of elimination in this complexity model. 

Oh 

< 

^ 1 Introduction 

(N 

Modern elimination theory starts with Kronecker's 1882 paper |Kro82] where the 
argumentation is essentially constructive, i.e., algorithmic. Questions of efficiency of 
algorithms become only indirectly and marginally addressed in this paper. However, 
c/5 later criticism of Kronecker's approach to algebraic geometry emphasized the algo- 

lHj rithmic inefficiency of his argumentation ( |Macl6j . |vdW50] ) . In a series of more 

^ recent contributions, that started with |CGH89] and ended up with |GHM + 98] . 

|GHH + 97j . |HMW01j and |GLS01] . it became apparent that this criticism is based 
on a too narrow interpretation of Kronecker's elimination method. In fact, these 
£C) contributions are, implicitly or explicitly, based on this method, nonwithstanding 

that they also contain other views and ideas coming from commutative algebra and 
algebraic complexity theory. 

A turning point was achieved by the combination of a new, global view of New- 
ton iteration with Kronecker's method ([G HM + 98j . |GHH + 97j ). The outcome was 
>- that elimination polynomials, although hard to represent by their coefficients, allow 

a reasonably efficient encoding by evaluation algorithms. This circumstance sug- 
gests to represent in elimination algorithms polynomials not by their coefficients 
but by arithmetic circuits (see |HS81] . |Kal88] and [FIK86] for the first steps in 
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this direction). This idea became fully realized by the "Kronecker" algorithm for 
the resolution of polynomial equation systems over algebraically closed fields. The 
algorithm was anticipated in |GHMP95j . |GHM+98j . jHMWOlj . [GLSOlj and imple- 
mented in a software package of identical name (see |Lecj ) . 

The results of the present paper imply that the complexity of the Kronecker 
algorithm is asymptotically optimal under reasonable architectural assumptions. In 
particular, we exhibit in Section 5J3 of this paper an infinite family of arithmetic 
input circuits encoding efficiently certain elimination problems such that any al- 
gorithm, designed by commonly accepted rules of software engineering, requires 
exponential time to solve these problems. In order to establish such a result we 
need a computation model. In this contribution we put main emphasis on this task. 

In the same spirit we prove in Section [4] that there exists an infinite family of 
parameter dependent elimination polynomials which, under a certain requirement 
of branching-freeness (called "robustness"), need arithmetic circuits of exponential 
size for their evaluation, whereas the circuit size of the corresponding elimination 
problems grows only polynomially. 

In the past, many attempts to show the non-polynomial character of the elimina- 
tion of just one existential quantifier block in the arithmetic circuit based elementary 
language over C, employed the reduction to the claim that an appropriate candidate 
family of specific polynomials was hard to evaluate (this approach was introduced 
in |HM93] and became adapted to the BSS model in |SS95j ). We give here the first 
example of such a family. 

We finish this paper with a discussion of the concept of an approximative algo- 
rithm and establish a link between this notion and our computation model. 



2 Concepts and tools from algebraic geometry 



In this section, we use freely standard notions and notations from commutative al- 
gebra and algebraic geometry. These can be found for example in |Lan93j . [ZS60J . 



|Kun85j and |Sha94] . In Sections 2.1 and 2.3 we introduce the notions and defini- 
tions which constitute our fundamental tool for the modelling of elimination prob- 
lems and algorithms. Most of these notions and their definitions are taken from 
jGHMSllj . 



2.1 Basic notions and notations 

For any n e N, we denote by A n := A n (C) the n-dimensional affine space C n 
equipped with its respective Zariski and Euclidean topologies over C. In algebraic 
geometry, the Euclidean topology of A n is also called the strong topology. We shall 
use this terminology only exceptionally. 

Let X±, . . . , X n be indeterminates over C and let X := (Xi, . . . , X n ). We denote 
by C[X] the ring of polynomials in the variables X with complex coefficients. 
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Let V be a closed affine subvariety of A n . As usual, we write dimV^ for the 
dimension of the variety V. 

For fi, . . . , f s , g G C[X] we shall use the notation {fi = 0, . . . , f s = 0} in order 
to denote the closed affine subvariety V of A n defined by fi, . . . , f s and the notation 
{/i = 0, . . . , f s = 0,g 7^ 0} in order to denote the Zariski open subset V g of V 
defined by the intersection of V with the complement of {g = 0}. Observe that V g 
is a locally closed affine subvariety of A n whose coordinate ring is the localization 

cm, of cm. 

We denote by I(V) : = {/ G C[X] : f(x) = for any x G V} the ideal of definition 
of V in C[X] and by C[V] := {if : V ->■ C ; there exists / G C[X] with <p(x) = 
f{x) for any x G V} its coordinate ring. Observe that C[V] is isomorphic to the 
quotient C-algebra C[V] := C[X]/I(V). If is irreducible, then C[V] is zero- 
divisor free and we denote by C(V) the field formed by the rational functions of V 
with maximal domain (C(V) is called the rational function field of V). Observe that 
C(V) is isomorphic to the fraction field of the integral domain C[V]. 

In the general situation where V is an arbitrary closed affine subvariety of A n , 
the notion of a rational function of V has also a precise meaning. The only point to 
underline is that the domain, say U, of a rational function of V has to be a maximal 
Zariski open and dense subset of V to which the given rational function can be 
extended. In particular, U has a nonempty intersection with any of the irreducible 
components of V. 

As in the case where V is irreducible, we denote by C(V) the C-algebra formed 
by the rational functions of V. In algebraic terms, C(V) is the total quotient ring 
of C[V] and is isomorphic to the direct product of the rational function fields of the 
irreducible components of V. 

Let be given a partial map <ft : V —■* W, where V and W are closed subvarieties 
of some affine spaces A n and A m , and let ... , 4> m be the components of <fi. With 
these notations we have the following definitions: 

Definition 1 (Polynomial map). The map <fi is called a morphism of affine varieties 
or just polynomial map if the complex valued functions <pi, . . . ,<f> m belong to C[V]. 
Thus, in particular, <p is a total map. 

Definition 2 (Rational map). We call <fi a rational map ofV to W , if the domain 
U of <f) is a Zariski open and dense subset of V and 0i, . . . , <p m are the restrictions 
of suitable rational functions of V to U. 

Observe that our definition of a rational map differs from the usual one in alge- 
braic geometry, since we do not require that the domain U of is maximal. Hence, in 
the case m :— 1, our concepts of rational function and rational map do not coincide 
(see also |GHMSll] l 
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2.2 Constructible sets and constructible maps 

Let Ai be a subset of some affine space A n and, for a given nonnegative integer m, 
let : Ai — ■* A m be a partial map. 

Definition 3 (Constructible set). We call the set Ai constructible if Ai is definable 
by a Boolean combination of polynomial equations. 

A basic fact we shall use in the sequel is that if Ai is constructible, then its 
Zariski closure is equal to its Euclidean closure (see, e.g., |Mum88] . Chapter I, §10, 
Corollary 1). In the same vein we have the following definition. 

Definition 4 (Constructible map). We call the partial map constructible if the 

graph of <p is constructible as a subset of the affine space A n x A m . 

We say that is polynomial if is the restriction of a morphism of affine varieties 
A n — > A m to a constructible subset Ai of A n and hence a total map from Ai to A m . 
Furthermore, we call a rational map of Ai if the domain U of is contained in 
Ai and is the restriction to Ai of a rational map of the Zariski closure Ai of Ai. 
In this case U is a Zariski open and dense subset of Ai. 

Since the elementary, i.e., first-order theory of algebraically closed fields with 
constants in C admits quantifier elimination, constructibility means just elementary 
definability. In particular, is constructible implies that the domain and the image 
of are constructible subsets of A n and A m , respectively. 

Remark 1. A partial map : Ai --• * A m is constructible if and only if it is piecewise 
rational. If is constructible there exists a Zariski open and dense subset U of Ai 
such that the restriction <fi\u of to U is a rational map of Ai (and of Ai ). 

For details we refer to [GHMSllJ, Lemma 1. 

2.3 Geometrically robust constructible maps 

The main mathematical tool of this paper is the notion of geometrical robustness 
which we are going to introduce now. 

Let Ai be a constructible subset of the affine space A n and let : Ai — > A m be 
a (total) constructible map with components 0i, . . . , m . 

We consider now the Zariski closure Ai of the constructible subset Ai of A n . 
Observe that is a closed affine subvariety of A n and that we may interpret C(A4) 
as a C[A4]-module (or C[A4] -algebra). 

Fix now an arbitrary point x of Ai. By %R X we denote the maximal ideal of 
coordinate functions of C[A4] which vanish at the point x. By CLVl]^ we denote 
the local C-algebra of the variety Ai at the point x, i.e., the localization of C[A4] at 
the maximal ideal %3l x . By ClAi)^ we denote the localization of the C[A^]-module 
C(M) at <m x . 
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Following Remark [TJ we may interpret (pi, ... , <f> m as rational functions of the 
affine variety Ai and therefore as elements of the total fraction ring C(Ai) of CLV1]. 
Thus C[jM][0i, . . . , 4> m ] and C[A^]gj^.[0i, . . . , 4> m ] are C-subalgebras of C(Ai) and 
C(Ai)^n x which contain C[.M] and CfA^]^, respectively. 

The following result establishes for constructible maps a bridge between a topo- 
logical and an algebraic notion. It will be fundamental in the context of this paper. 

Theorem— Definition 2. Let notations and assumptions be as before. We call the 
constructible map <fi : Ai — > A m geometrically robust if cf) is continuous with respect to 
the Euclidean topologies of Ai and A m or equivalently, if <pi, . . . , <p m , interpreted as 
rational functions of the affine variety Ai, satisfy at any point x G Ai the following 
two conditions: 

(i) £\Ai]m x [fa, ■ ■ ■ ></>m] is a finite C[Ai]<xn x -module. 

(ii) C[Ai]tm x [fa, ■ ■ ■ , 4>m} is a local C[Ai](xn x - algebra whose maximal ideal is gener- 
ated by M x and fa - fa(x), ...,<p m - <j) m (x). 

For a proof of this result, which is based on Zariski's Main Theorem ( |Ive73j . 
§IV.2) we refer to jHKRllj (see also |CGH+03| and [GHMS11] ). 

From the topological definition of a geometrically robust constructible map one 
deduces immediately the following statement. 

Corollary 3. If we restrict a geometrically robust constructible map to a con- 
structible subset of its domain of definition we obtain again a geometrically robust 
map. Moreover the composition and the cartesian product of two geometrically robust 
constructible maps are geometrically robust. The geometrically robust constructible 
functions form a commutative C-algebra which contains the polynomial functions. 

The origin of the concept of a geometrically robust map can be found, implic- 
itly, in [GH01J. It was introduced explicitly for constructible maps with irreducible 
domains of definition in |GHMSlT] , where it is used to analyze the complexity char- 
acter of multivariate Hermite-Lagrange interpolation. 

3 Robust parameterized arithmetic circuits 

We shall use freely standard concepts from algebraic complexity theory which can 
be found in |BCS97j . 

Let us fix natural numbers n and r, indeterminates X±, . . . , X n and a non-empty 
constructible subset Ai of A r . By tti, . . . , 7T r we denote the restrictions to Ai of the 
canonical projections A r — > A 1 . 

A (by Ai ) parameterized arithmetic circuit f3 (with basic parameters m, . . . ,ir r 
and inputs X±, . . . , X n ) is a labelled directed acyclic graph (labelled DAG) satisfying 
the following conditions: 

each node of indegree zero is labelled by a scalar from C, a basic parameter m, . . . ,ir r 
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or a input variable Xi, . . . , X n . Following the case, we shall refer to the scalar, basic 
parameter and (standard) input nodes of (3. All other nodes of (3 have indegree 
two and are called internal. They are labelled by arithmetic operations (addition, 
subtraction, multiplication, division). A parameter node of (3 depends only on scalar 
and basic parameter nodes, but not on any input node of (3 (here "dependence" refers 
to the existence of a connecting path). An addition or multiplication node whose 
two ingoing edges depend on an input is called essential. The same terminology is 
applied to division nodes whose second argument depends on an input. Moreover, 
at least one circuit node becomes labelled as output. Without loss of generality we 
may suppose that all nodes of outdegree zero are outputs of (3. 

We consider (3 as a syntactical object which we wish to equip with a certain 
semantics. In principle there exists a canonical evaluation procedure of (3 assigning 
to each node a rational function of Ai x A" which, in case of a parameter node, may 
also be interpreted as a rational function of Ai. In either situation we call such a 
rational function an intermediate result of (3. 

The evaluation procedure may fail if we divide at some node an intermediate 
result by another one which vanishes on a Zariski dense subset of a whole irreducible 
component of M x A". If this occurs, we call the labelled DAG (3 inconsistent, 
otherwise consistent. 

If nothing else is said, we shall from now on assume that (3 is a consistent 
parameterized arithmetic circuit. The intermediate results associated with output 
nodes will be called final results of (3. 

We call an intermediate result associated with a parameter node a parameter of 
(3 and interpret it generally as a rational function of Ai. A parameter associated 
with a node which has an outgoing edge into a node which depends on some input 
of (3 is called essential. In the sequel we shall refer to the constructible set Ai as 
the parameter domain of (3. 

We consider (3 as a syntactic object which represents the final results of (3, i.e., 
the rational functions of Ai x A ra assigned to its output nodes. 

Now we suppose that the parameterized arithmetic circuit (3 has been equipped 
with an additional structure, linked to the semantics of (3. We assume that for 
each node p of (3 there is given a total constructible map Ai x A™ — >■ A 1 which 
extends the intermediate result associated with p. In this way, if (3 has K nodes, we 
obtain a total constructible map Q : Ai x A n — > A K which extends the rational map 
Ai x A™ --■» A^ given by the labels at the indegree zero nodes and the intermediate 
results of (3. 

Definition 5 (Robust circuit). Let notations and assumptions be as before. The 
pair ((3, Q) is called a robust parameterized arithmetic circuit if the constructible 
map Q is geometrically robust. 

Observe that the above rational map Ai x A™ --■» A K can be extended to at 
most one geometrically robust constructible map 0:^x4"-^ A K . Therefore we 
shall apply from now on the term "robust" also to the circuit (3. 
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Robust parameterized arithmetic circuits may be restricted as follows: 
Let TV be a constructible subset of M. and suppose that ((3, Q) is robust. Then 
Corollary [3] implies that the restriction £l\fj X A n of the constructible map to J\fx A n 
is still a geometrically robust constructible map. 

This implies that (/3, Q) induces a by M parameterized arithmetical circuit (3jy. 
Observe that /3/v may become inconsistent. If /3/v is consistent then ((3j*j, fl\jVxA n ) is 
robust. The nodes where the evaluation of /3/v fails correspond to divisions of zero 
by zero which may be replaced by so called approximative algorithms having unique 



limits (see |HKRllj . Section 3.3.2 and Section 5.6 below). These limits are given by 
the map fl\j^ X A n - We call (f3j^, Q\jVxA n ), or simply (3^, the restriction of (/3, O) or (3 
to M. 

We say that the parameterized arithmetic circuit /3 is totally division-free if any 
division node of /3 corresponds to a division by a non-zero complex scalar. 

We call essentially division-free if only parameter nodes are labelled by divi- 
sions. Thus the property of (3 being totally division-free implies that /3 is essentially 
division-free, but not vice versa. Moreover, if /3 is totally division-free, the rational 
map given by the intermediate results of (3 is polynomial and therefore a geometri- 
cally robust constructible map. Thus, any by Ai parameterized, totally division-free 
circuit is in a natural way robust. 

In the sequel, we shall need the following elementary fact. 

Lemma 4. Let notations and assumptions be as before and suppose that the pa- 
rameterized arithmetic circuit (3 is robust. Then all intermediate results of (3 are 
polynomials in X 1 , . . . ,X n over the C-algebra of geometrically robust constructible 
functions defined on A4 . 



For a proof of Lemma [4] we refer to |HKRllj . Section 3.1. 

The statement of this lemma should not lead to confusions with the notion of 
an essentially division-free parameterized circuit. We say just that the intermediate 
results of f3 are polynomials in X±, . . . , X n and do not restrict the type of arithmetic 
operations contained in (3 (as we did defining the notion of an essentially division- 
free parameterized circuit). 

To our parameterized arithmetic circuit f3 we may associate different complexity 
measures and models. In this paper we shall mainly be concerned with sequential 
computing time, measured by the size of (3. Here we refer with "size" to the number 
of internal nodes of (3 which count for the given complexity measure. Our basic 
complexity measure is the non-scalar one (also called Ostrowski measure) over the 
ground field C. This means that we count, at unit costs, only essential multiplica- 
tions and divisions (involving basic parameters or input variables in both arguments 
in the case of a multiplication and in the second argument in the case of a division), 
whereas C-linear operations are free (sec [BCS97J for details). 
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3.1 Operations with robust parameterized arithmetic cir- 
cuits 

3.1.1 The operation join 

Let 71 and 72 be two robust parameterized arithmetic circuits with parameter do- 
main Ai and suppose that there is given a one-to-one correspondence A which 
identifies the output nodes of 71 with the input nodes of 72 (thus they must have 
the same number). Using this identification we may now join the circuit 71 with the 
circuit 72 in order to obtain a new parameterized arithmetic circuit 72 *a 7i with 
parameter domain Ai. The circuit 72 *a 7i has the same input nodes as 71 and the 
same output nodes as 72 and one deduces easily from Lemma [4] and Corollary [3] 
that the circuit 72 *a 7i is robust and represents a composition of the rational maps 
defined by 71 and 72, if 72 *a 7i is consistent. The (consistent) circuit 72 *a 7i is 
called the (consistent) join of 71 with 72. 

Observe that the final results of a given robust parameterized arithmetic circuit 
may constitute a vector of parameters. The join of such a circuit with another robust 
parameterized arithmetic circuit, say (3, is again a robust parameterized arithmetic 
circuit which is called an evaluation of (3. Hence, mutatis mutandis, the notion of 
join of two routines includes also the case of circuit evaluation. 

3.1.2 The operations reduction and broadcasting 

We describe now how, based on its semantics, a given parameterized arithmetic 
circuit (3 with parameter domain Ai may be rewritten as a new circuit over Ai 
which computes the same final results as (3. 

The resulting two rewriting procedures, called reduction and broadcasting, will 
neither be unique, nor generally confluent. To help understanding, the reader may 
suppose that there is given an (efficient) algorithm which allows identity checking 
between intermediate results of (3. However, we shall not make explicit reference to 
this assumption. We are now going to explain the first rewriting procedure. 

Suppose that the parameterized arithmetic circuit (3 computes at two different 
nodes, say p and p', the same intermediate result. Assume first that p neither 
depends on p', nor p' on p. Then we may erase p' and its two ingoing edges (if p' is 
an internal node) and draw an outgoing edge from p to any other node of (3 which is 
reached by an outgoing edge of p'. If p' is an output node, we label p also as output 
node. Observe that in this manner a possible indexing of the output nodes of (3 may 
become changed but not the final results of (3 themselves. 

Suppose now that p' depends on p. Since the DAG (3 is acyclic, p does not 
depend on p' . We may now proceed in the same way as before, erasing the node p'. 

Let j3' be the parameterized arithmetic circuit obtained, as described before, by 
erasing the node p'. Then we call (3' a reduction of (3 and call the way we obtained f3' 
from (3 a reduction step. A reduction procedure is a sequence of successive reduction 
steps. 
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One sees now easily that a reduction procedure applied to produces a new- 
parameterized arithmetic circuit (3* (also called a reduction of 0) with the same basic 
parameter and input nodes, which computes the same final results as (3 (although 
their possible indexing may be changed). Moreover, if is a robust parameterized 
circuit, then (3* is robust too. Observe also that in the case of robust parameterized 
circuits our reduction commutes with restriction. 

Now we introduce the second rewriting procedure. 

Let assumptions and notations be as before and let be given a set P of nodes 
of f3 and a robust parameterized arithmetic circuit 7 with parameter domain Ai 
and #P input nodes, namely for each p G P one which becomes labelled by a new 
input variable Y p . We obtain a new parameterized arithmetic circuit, denoted by 
7 *p (3, when we join 7 with f3, replacing for each p G P the input node of 7, which 
is labelled by the variable Y p , by the node p of (3. The output nodes of (3 constitute 
also the output nodes of 7 *p (3. Thus (3 and 7 *p (3 compute the same final results. 
Observe that 7 *p (3 is robust if it is consistent. We call the circuit 7 *p (3 and all its 
reductions broadcastings of (3. Thus broadcasting a robust parameterized arithmetic 
circuit means rewriting it using only valid polynomial identities. 

If we consider arithmetic circuits as computer programs, then reduction and 
broadcasting represent a kind of program transformations. 

4 A family of hard elimination polynomials 

As a major result of this paper we are now going to exhibit an infinite family of pa- 
rameter dependent elimination polynomials which require essentially division-free, 
robust parameterized arithmetic circuits of exponential size for their evaluation, 
whereas the circuit size of the corresponding input problems grows only polyno- 
mially. The proof of this result, which is absolutely new in his kind, is astonishly 
elementary and simple. 

Let T, Ui, . . . , U n and X 1; . . . , X n be indeterminates and let U := (U%, . . . , U n ) 
and X : = (X 1; . . . , X n ). Consider for given n6N the polynomial : = ^ 1<i<n 2 l ~ 
Xi + T"r[ 1<i<n (l + {Ui — l)Xi). Observe that can be evaluated using n — 1 
non-scalar multiplications involving X%, . . . , X n . 

The set O := {Ei<,<„ 2^ + t Ui^ni 1 + («< - W; (*, "!>■■■> "n) e A« +1 } 
is contained in a finite-dimensional C-linear subspace of C[X] and therefore O and 
its closure O are constructible sets. 

From |GHMSll] . Section 3.3.3 we deduce the following facts: 
there exist K := 16n 2 + 2 integer points £1, e Z n of bit length at most 4n 

such that for any two polynomials f,g £ O the equalities f((,k) = g(£,k), 1 < k < K , 
imply f = g. Thus the polynomial map H : O — > A defined for / G O by 
S(/) := (/(£i), is injective. Moreover M. := S(O) is an irreducible 
constructible subset of and we have Ai = Finally, the constructible map 

:= which maps A4 onto (9 and Ai onto (9, is a restriction of a geometrically 
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robust map and therefore by Corollary [3] itself geometrically robust. 

For e G {0, l} n we denote by 4> £ the map M. — > A 1 which assigns to each point v G 
A4 the value 0(f) (e). From Corollary [3] we conclude that e is a geometrically robust 
constructible function which belongs to the function field C(Ai) of the irreducible 
algebraic variety Ai. 

Observe that for t G A 1 and u G A n the identities (p t (E(H^(t, u, X))) = 
(j)(E(H^(t,u,X)))(e) = ((H- 1 oH)(ffW(t,ii,I)))(e) = #»(t, u, e) hold. 

Let p( n ) := riee{o i} n 0^ ~~ ^«)- Then p( n ) is a geometrically robust constructible 
function which maps M x A 1 (and hence MxA 1 ) into A 1 . 

Consider now the polynomial := ]\ eem}n (Y - (T ,U , e)) = n o <j<2"-i( y - 

i<i<n ^* '*))> where [j]i denotes the i-th digit of the binary representation 
of the integer j, < j < 2 n — 1, 1 < i < n. We have for t G A 1 and 116 A" the 
identities 

pM(E(HV>(t, U ,X)),Y)= (Y-4> e (E(H^\t,u,X)))) = 

J] (F-ifW(t,n,6)) = FM(t, M ,r) 

e£{0,l}" 

Let Si, . . . , Sk be new indeterminates and observe that the existential first order 
formula of the elementary theory of C, namely 

(3Xi) . . . (3X n )(3T)(3fA) . . . (3f/ n )(X 1 2 - X 1 = A ■ • • A X 2 n - X n = OA 
/\ Sj = H^\T,U,Q AY = H {n) (T,U,X)) (2) 

l<j<K 

describes the constructible subset {(s,y) G A K+1 ; s G A4,y G A 1 , p( n \s, y) = 
0} of A K+1 . Moreover, P^ is the greatest common divisor in C(.M)[F] of all 
polynomials of C[.M][F] which vanish identically on the constructible subset of A^ +1 
defined by the formula ([2]). Hence P^ G C(.M)[y] is a (parameterized) elimination 
polynomial. 

Observe that the polynomials contained in the formula ^ can be represented 
by a totally division-free arithmetic circuit of size 0(n 3 ). Therefore, the formula 
§2§ is also of size 0(n 3 ). 

Theorem 5. Let notations and assumptions be as before and let 7 be an essentially 
division- free, robust parameterized arithmetic circuit with domain of definition M. 
such that 7 evaluates the elimination polynomial P( n \ 

Then 7 performs at least f2(2a ) essential multiplications and at least Q(2 n ) mul- 
tiplications with parameters. 

Proof. We fix the natural number n. Let us write H := = ^2 1<i<n T~ 1 X i + 

Y\i <i<n T{Ui — l)Xi as a polynomial in the main indeterminates Xi, . . . ,X n with 
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coefficients 9 Klr .. )Kn G C[T, U], n\, . . . , n n G {0, 1}, namely 

H= Yl 9 Kl ,... >Kn X^,...,X^. 
Ki,...,K n e{o,i} 

Observe that for Ki, . . . , K n G {0, 1} the polynomial 5 K1) ... lKn (0, U) G C[Z7] is of degree 
at most zero, i.e., a constant complex number, independent of U\, . . . , U n . 

Let 9 := (^«i,.„,K„)/ci,...,/« n e{Q ) i} an d observe that the vector 0(0, U) is a fixed point 
of the affine space A 2 ™. We denote by SOT the vanishing ideal of the C-algebra C[9] at 
this point. We interpret 9 as a geometrically robust constructible map A rt+1 — > A 2 " 
with (constructible) image T. 

Let us write F := = rio<j<2"-i( r ~ C? + T Ili<i<n ^O) as a polynomial 
in the main indeterminate Y with coefficients y? K G C[T, £/], 1 < k < 2 n , namely 
F = Y 2n + i Pl Y 2n - 1 + --- + cp 2n . 

Observe that for 1 < k < 2 n the polynomial (p K (0, U) G C[U] is of degree at most 
zero. Let X K := (p K (0,U), X := (A K )i< K < 2 « and := (<^ k )i<k<2™- Observe that A is a 
fixed point of the affine space A 2 ™. 

Let v : A™ +1 — > A K be the polynomial map defined for t G A 1 and u G A n by 
is(t, u) := E(H(t, u, X)) = (H(t, u, £i), . . . , H(t, u, £,k))- Observe that there exists a 
geometrically robust constructible map o : T — > A K such that ao9 = v holds. Since 
by assumption the parameterized arithmetic circuit 7 is essentially division-free and 
robust, there exists a geometrically robust constructible map ip defined on M. such 
that the entries of ip constitute the essential parameters of the circuit 7. Moreover, 
for m being the number of components of i/j, there exists a vector u of m-variate 
polynomials over C such that the entries of u)(ip) — ou 01/j become the coefficients of 
the elimination polynomial P := = rLe{o i} n 0^ ~ 

One sees easily that there exists a totally division-free ordinary arithmetic circuit 
7' which evaluates the polynomials H(T, U, £1), . . . , H(T, U, 

The join 7 * 7' of 7' with 7 at the basic parameter nodes of 7 is an essentially 
division-free robust parameterized circuit with domain of definition A™ +1 which by 
([!]) evaluates the polynomial F(T } U, Y) := F^ n \T } U, Y). The entries of the vector 
v := if; o v constitute the essential parameters of the circuit 7*7' and the entries of 
ujou = u)oil)ou become by ([!]) the coefficients of the polynomial F(T, U, Y) with 
respect to Y . So we have <p — u o v. 

Taking into account v = ipov = ipoao9, Theorem-Definition|2](z) and -01 l.MSTTj . 
Corollary 12 we conclude that the entries of v are polynomials of C[T, U] which are 
integral over the local C-subalgebra C[0]an of C(T, U). 

Let G C[T, U] be such an entry. Then there exists an integer s and polynomials 
Q-O; Oil , . . . , CL S G C[9] with do <£Wl such that the algebraic dependence relation 

a n s + ai/r* -1 H + a s = (3) 

is satisfied in C[T, U]. From ^ we deduce the algebraic dependence relation 

a (0, U)n(0, U) s + ai(0, U)n{0, U) 8 ' 1 + ■■■ + a s (0, U) = (4) 
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in C[U]. 

Since the polynomials do, OX) ■ • • , CLs 

belong to C[9] and 9(0, U) is a fixed point of 
A 2 ", we conclude that a := a (0, U), a\ := ai(0, U), . . . , a s := a s (0, U) are complex 
numbers. Moreover, a VJt implies a Q ^ 0. 

Thus Q may be rewritten into the algebraic dependence relation 

aoMO, U) s + ai/i(0, f/) 3 " 1 + • • • + a s = (5) 

in C[U] with a ^ 0. 

This implies that the polynomial fi(0,U) of C[U] is of degree at most zero. 
Therefore w := z>(0, [/) is a fixed point of the affine space A m . 

Recall that A = (A K )i< K <2" with X K := <^ K (0, U), 1 < k < 2 n , is a fixed point of 
the affine space A 2 ". 

From |CGH + 03j . Lemma 6 we deduce that for 1 < n < 2 n the coefficient cp K of 
F is an element of C[T, U] of the form 

<Pk = A K + TL K + terms of higher degree in T (6) 

where L\, . . . , L 2 n G C[U] are C-linearly independent. 

Consider now an arbitrary point u G A™ and let e u : A 1 — > A m and S u : A 1 — > A 2 " 
be the polynomial maps defined for t G A 1 by e„(t) := v{t,u) and := ip(t,u). 

Then we have e u (0) = z>(0,m) = u> and 5 M (0) = f(0,u) = A, independently of u. 
Moreover, from ip = oj o u we deduce (5 U = u o e u . 

Thus ([6]) implies 

. . . ,L 2 n(u)) = ^(0,«) = ^(0) = (^) w «(0)), (7) 

where (Du) w denotes the (first) derivative of the m-variate polynomial map u 
at the point w G A' m and S' u (0) and e' u (0) are the derivatives of the parameterized 
curves 5 U and e u at the point G A 1 . We rewrite now ([7]) in matrix form, replacing 
(Du) w by the corresponding transposed Jacobi matrix M G A mx2 " and S' u (0) and 
e' u (0) by the corresponding points of A 2 ™ and A m , respectively. 

Then ^ takes the form 

(L 1 (u),...,L 2n (u)) = e' u (0)M, (8) 

where the complex (m x 2™)-matrix M is independent of u. 

Since the polynomials L±, . . . ,L 2 ™ G C[U] are C-linearly independent, we may 
choose points Ui, . . . G A™ such that the complex (2 n x 2 n )-matrix 

N := (L K (M i ))i< ijK < 2 n 

has rank 2™. 

Let K be the complex (2™ x m)-matrix whose rows are 6^(0), . . . , e' u n (0). 
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Then rt8j) implies the matrix identity 



N = K ■ M. 

Since N has rank 2™, the rank of the complex (m x 2 ra )-matrix M is at least 2 n . 
This implies 

m > 2 n . (9) 

Therefore the circuit 7 contains m > 2 n essential parameters. 

Let L be the number of essential multiplications executed by the parameterized 
arithmetic circuit 7 and let V be the total number of multiplications of 7, excepting 
those by scalars from C. Then, after a well-known standard rearrangement [PS73J 
of 7, we may suppose without loss of generality, that there exists a constant c > 
(independent of the input circuit 7) such that L > cmz and L' > cm holds. 

Let L be the number of essential multiplications executed by the parameterized 
arithmetic circuit 7 and let V be the total number of multiplications of 7, excepting 
those by scalars from C. Then, after a well-known standard rearrangement |PS73] 
of 7, we may suppose without loss of generality, that there exists a constant c > 
(independent of the input circuit 7 and the procedure ^4) such that L > cm 2 and 
V > cm holds. 

From the estimation ^ we deduce now that the circuit 7 performs at least 
Q(2?) essential multiplications and at least Q(2 n ) multiplications, including also 
multiplications with parameters. This finishes the proof of the theorem. □ 

Theorem [5] is essentially contained in the arguments of the proof of |GH01j . 
Theorem 5 and |CGH + 03j . Theorem 4. 

Observe that a quantifier-free description of Ai by means of circuit represented 
polynomials, together with an essentially division-free, robust parameterized arith- 
metic circuit 7 with domain of definition Ai, which evaluates the elimination polyno- 
mial P^ captures the intuitive meaning of an algorithmic solution of the elimination 
problem described by the formula (|2|, when we restrict our attention to solutions 
of this kind and minimize the number of equations and branchings. In particular 
the circuit 7 can be evaluated for any input point (s, y) with s G Ai and y G C 
and the intermediate results of 7 are polynomials of C(.M)[Y] whose coefficients are 
geometrically robust constructible functions defined on Ai. 

With respect to the indeterminate Y, the coefficients of the polynomial G 
C(.M)[F] are geometrically robust constructible functions of the parameter domain 
Ai. In order to consider as an elimination polynomial as we did, the reader 
might expect that the coefficients of P^ should belong, for any point s G Ai, to 
the local ring of Ai at s. This would be true if the algebraic variety Ai would be 
normal at any s G Ai (see jGHMSllj . Corollary 12). From |CGH+03j . Corollary 3 
we deduce that the variety Ai is definitely not normal. This leads us to the question 
how elimination polynomials should look like when the closure of the parameter 
domain is not normal. 
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In order to elucidate this question we shall consider the following general situa- 
tion. It turns out that the requirement that the coefficients of elimination polyno- 
mials should be geometrically robust constructible functions is quite natural. 

Let ip : V — > W be a finite surjective morphism of irreducible affine varieties V 
and W over C such that there exists a coordinate function y G C[V] with C[V] = 
C[W][y]. Let d := [C(V) : C(W / )] be the degree of <p and suppose that for any 
point w G W the cardinality of the fiber (f~ 1 (w) is exactly d. Finally, let Y be a 
new indeterminate and F := Y d + <p d „ 1 Y d - 1 + --- + (p G C(W)[Y] the minimal 
polynomial of y. Observe that coefficients of F, namely (po, . . . , (fd-i G C(W), are 
integral over C[W]. We are now going to discuss a condition under which F may be 
considered as an elimination polynomial. This condition will imply that <po, . . . , (fid-i 
are geometrically robust constructible functions. 

We shall use the following abbreviations: 
A := C[W], B := A[<p , . . . , C := A[y], D := B[y\. We have the following 

commutative diagram of integral C-algebra extensions: 



B 




A D 




C 



Observe that D is isomorphic to B\Y]/ b[y]-f and in particular a free 5-module of 
rank d. 

Proposition 6. Suppose that for any maximal ideal m of A the canonical C-algebra 
homomorphism A/ m — >■ C/ m c is unramified (\lve.73j , Chapter I) and that m is con- 
tained in at most d maximal ideals of D (thus, intuitively, F is an elimination 
polynomial) . Then (po, . . . , (fa-i are geometrically robust constructible functions of 
W. 

Proof. Let m be an arbitrary maximal ideal of A. Since A — > B is an integral ring 
extension, we deduce from Theorem-Definition [2] that it suffices to show that there 
exists a single maximal ideal of B which contains m. Our assumptions yield a 
commutative diagram 

A/ m C/ m c 

i Hi 

c " c d 
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Taking into account C = A[y] we conclude that there exists a monic polynomial 
G G A[y] of degree d with discriminant p G A such that C/ m c is isomorphic to A[Y] 
divided by the ideal generated m and G and such that p does not belong to m. 

Let 91 be an arbitrary maximal ideal of B which contains m and let F and G 
be the images of F and G in B/<ji[Y]. Then we have DJ<rd — B/ra[Y]/ B ,r Y -i.p and 
therefore F divides G in B/sr\Y\. From d = degF = degG and the fact that F and 
G are monic we deduce F = G. Since the discriminant p of G does not belong to 
m we have p ^ 91 and therefore the polynomial F is separable. Thus we obtain a 
commutative diagram 

B/m *■ D/<xd 

i i 

C C d 

and in particular the canonical C-algebra homomorphism B/m — » -D/oto is un- 
ramified. Hence the number of maximal ideals of D which contain 91 is exactly d. 
By assumption there are at most d maximal ideals of D containing m. Therefore any 
such ideal must contain 91. Since B — > D is an integral ring extension, we conclude 
that 91 is the unique maximal ideal of B which contains m. □ 

5 A computation model with robust parameter- 
ized arithmetic circuits 

This section is devoted to a deeper understanding of the assumptions which lead to 
Theorem [5j To this end we introduce a computation model which will be comprehen- 
sive enough to capture the essence of all known circuit based elimination algorithms 
in effective algebraic geometry and, mutatis mutandis, also of all other (linear alge- 
bra and truncated rewriting) elimination procedures (see |Mor03] . |Mor05j . and the 
references cited therein, and for truncated rewriting methods especially [DFGS91]). 

The elimination problem and polynomial of Section [4] were somewhat artificial. 
We shall show that the conclusions of Theorem [5] are still valid for much more natural 
elimination problems and polynomials if we restrict the notion of algorithm to the 
computation model we are going to introduce in this section. 

The routines of our computation model will transform a given robust parameter- 
ized arithmetic (input) circuit into another robust parameterized arithmetic (output) 
circuit such that both circuits have the same parameter domain. 

In the sequel, we shall use ordinary arithmetic circuits over C as generic com- 
putations [BCS97J (also called computation schemes in |Hei89] ) . The indegree zero 
nodes of these arithmetic circuits are labelled by scalars and parameter and input 
variables. 

The aim is to represent different parameterized arithmetic circuits of similar size 
and appearance by different specializations (i.e., instantiations) of the parameter 
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variables in one and the same generic computation. For a suitable specialization 
of the parameter variables, the original parameterized arithmetic circuit may then 
be recovered by an appropriate reduction process applied to the specialized generic 
computation. 

This alternative view of parameterized arithmetic circuits will be fundamental 
for the design of routines of the computation model we are going to describe in 
this section. The routines of our computation model will operate on robust param- 
eterized arithmetic circuits and their basic ingredients will be subroutines which 
calculate parameter instances of suitable, by the model previously fixed, generic 
computations. These generic computations will be organized in finitely many fam- 
ilies which will only depend on a constant number of discrete parameters. These 
discrete families constitute the basic building block of our computation model (see 
|HKR11] . Section 3.2 for details about generic computation). 

In the sequel we shall distinguish sharply between the notions of input variable 
and parameter and the corresponding categories of circuit nodes. 

Input variables, called "standard", will occur in parameterized arithmetic circuits 
and generic computations. The input variables of generic computations will appear 
subdivided in three sorts, namely as "parameter" , "argument" and "standard" input 
variables. 

The computation model we are going to introduce now will assume different 
shapes, each shape being determined by a finite number of a priori given discrete 
(i.e., by tuples of natural numbers indexed) families of generic computations. The 
labels of the inputs of the ordinary arithmetic circuits which represent these generic 
computations will become subdivided into parameter, argument and standard input 
variables. We shall use the letters like U, U', U", . . . and W, W, W" to denote vec- 
tors of parameters, Y,Y' ,Y" , . . . and Z, Z' , Z" to denote vectors of argument and 
X, X', X", ... to denote vectors of standard input variables. 

We shall not write down explicitly the indexings of our generic computations by 
tuples of natural numbers. Generic computations will simply be distinguished by 
subscripts and superscripts, if necessary. 

Ordinary arithmetic circuits of the form 

Rx^W^XW), R X2 (W 2 ;X^), ... 
R' Xl (Wy,XW), R' X2 (W 2 ,;X^), ... 

represent a first type of a discrete family of generic computations (for each variable 
Xi,X 2 , . . . ,X n we suppose to have at least one generic computation). Other types 
of families of generic computations are of the form 

R + (W;U,Y;X), 
R.,(W',U,Y',X), 
R add {W;Y,Z;X), 
R mult (W;Y,Z;X), 



R' + {W';U',Y';X'), 
R[ f (W';U',Y';X'), 
R' add (W';Y',Z';X'), 
R' mult (W';Y',Z';X'), 



R\(W";U",Y";X") 
R" / (W";U",Y";X") 
K dd iW";Y",Z";X") 
R'Lu(W";Y",Z";X") 
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and 



R dlv (W;Y,Z;X), R' div (W';Y>,Z';X'), R'^ V (W"; Y", Z"; X") .... 

Here the subscripts refer to addition of, and multiplication or division by a 
parameter (or scalar) and to essential addition, multiplication and division. A final 
type of families of generic computations is of the form 

R(W;Y;X), R'(W';Y';X'), R"(W"; Y"\ X"), . . . 

The inputs of the circuits handled by our computation model will only consist of 
standard variables. 

From now on we have in mind a previously fixed shape when we refer to the 
computation model we are introducing. We start with a given finite set of discrete 
families of generic computations which constitute a shape as described before. 

5.1 The notions of well behavedness under restrictions, 

isoparametricity and well behavedness under reductions 

A fundamental issue is how we recursively transform a given input circuit into an- 
other one with the same parameter domain. During such a transformation we make 
an iterative use of previously fixed generic computations. On their turn these deter- 
mine the corresponding recursive routine of our branching-free computation model. 

We consider again our input circuit /3. We suppose that we have already chosen 
for each node p, which depends at least on one of the input variables X±, . . . , X n , a 
generic computation 

R%(W P ;X^), 

R ( f(W p ,U p ,Y p ,X^), 
R^(W p ,U p ,Y p ,X^), 

R { a P l(W P ;Y p ,Z p ;X^), 

r ( :Uw p ,y p ,z p ,x^), 

R^ V {W P -Y P ,Z P -X^), 

and that this choice was made according to the label of p, namely X iy 1 < i < n, 
or addition of, or multiplication or division by an essential parameter, or essential 
addition, multiplication or division. Here we suppose that U p is a single variable, 
whereas W p , Y p , Z p and X^' may be arbitrary vectors of variables. 

Furthermore, we suppose that we have already precomputed for each node p 
of (3, which depends at least on one input, a vector w p of geometrically robust 
constructible functions defined on Ai. If p is an input node we assume that w p is a 
vector of complex numbers. Moreover, we assume that the length of w p equals the 
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length of the variable vector W p . We call the entries of w p the parameters at the 
node p of the routine A applied to the input circuit j3. 

We are now going to develop the routine A step by step. The routine A takes 
over all computations of (3 which involve only parameter nodes, without modifying 
them. 

Consider an arbitrary internal node p of (3 which depends at least on one input. 
The node p has two ingoing edges which come from two other nodes of f3, say p\ 
and p 2 . Suppose that the routine A, on input (3, has already computed two results, 
namely F pi and F p2 , corresponding to the nodes pi and p 2 - Suppose inductively that 
these results are vectors of polynomials depending on those standard input variables 
that occur in the vectors of the form X^ p \ where p' is any predecessor node of p. 
Furthermore, we assume that the coefficients of these polynomials constitute the 
entries of a geometrically robust, constructible map defined on M.. Finally we 
suppose that the lengths of the vectors F pi and Y p (or U p ) and F p2 and Z p coincide. 

The parameter vector w p of the routine A forms a geometrically robust, con- 
structible map defined on A4, whose image we denote by K p . Observe that JC P is a 
constructible subset of the affine space of the same dimension as the length of the 
vectors w p and W p . Denote by n p the vector of the restrictions to K p of the canon- 
ical projections of this affine space. We consider /C p as a new parameter domain 
with basic parameters n p . For the sake of simplicity we suppose that the node p is 
labelled by an essential multiplication. Thus the corresponding generic computation 
has the form 

R ( ±t(W P ,Y p ,Z p ,X^). (10) 
Let the specialized generic computation 

be the by /C„ parameterized arithmetic circuits obtained by substituting in the 



generic computation (10) for the vector of parameter variables W p the basic pa- 
rameters hi p . At the node p we shall now make the following requirement on the 
routine A applied to the input circuit (3: 

(A) The by K, p parameterized arithmetic circuit of Rl^ ult (hi p ;Y p , Z p ; X^), should 
be consistent and robust. 

Observe that the requirement (A) is automatically satisfied if all the generic 
computations of our shape are realized by totally division-free ordinary arithmetic 
circuits. 

Assume now that the routine A applied to the circuit (3 satisfies the requirement 
(A) at the node p of (3. 

Recall that we assumed that the node p is labelled by an essential multiplication 
and that the vectors F Pl and Y p and F P2 and Z p have the same length. Joining with 
the generic computation 

R ( ±u(W p ;Y p ,Z p ,X^) 
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at W p , Y p and Z p the previous computations of w p , F Pl and F P2 we obtain a param- 
eterized arithmetic circuit with parameter domain Ai, whose final results are the 
entries of a vector which we denote by F p . 

One deduces easily from our assumptions on w p , F pi and F p2 and from the re- 
quirement (A) in combination with Lemma [4] and Corollary [3j that the resulting 
parameterized arithmetic circuit is robust if it is consistent. The other possible la- 
bellings of the node p by arithmetic operations are treated similarly In particular, in 
case that p is an input node labelled by the variable X il 1 < i < n, the requirement 
(A) implies that the ordinary arithmetic circuit R^\w p ;X^) is consistent and ro- 
bust and that all its intermediate results are polynomials in over C (although 
Rx\w p ; X^) may contain divisions). 

We call the recursive routine A (on input (3) well behaved under restrictions if the 
requirement (A) is satisfied at any node p of (3 which depends at least on one input 
and if joining the corresponding generic computation with w p , F pi and F P2 produces 
a consistent circuit (observe that this last condition is automatically satisfied when 
the specialized generic computation of (A) is essentially division-free). If the routine 
A is well behaved under restrictions, then A transforms step by step the input circuit 
(3 into another consistent robust arithmetic circuit, namely *4(/3), with parameter 
domain A4. 

As a consequence of the recursive structure of A((3), each node p of generates a 
subcircuit of A((3) which we call the component of A((3) generated by p. The output 
nodes of each component of A((3) form the hypernodes of a hypergraph Ha(/3) whose 
hyperedges are given by the paths connecting the nodes of A((3) contained in distinct 
hypernodes of 'H.j^py The hypergraph Ha(P) m8b Y be shrunk to the DAG structure 
of (3 and therefore we denote the hypernodes of %A{p) i n the same way as the nodes 
of j3. Notice that well behavedness under restrictions is in fact a property which 
concerns the hypergraph %AW)- 

We call A a (recursive) parameter routine if A does not introduce new standard 
variables. In the previous recursive construction of the routine A, the parameters 
at the nodes of (3, used for the realization of the circuit A((3), are supposed to be 
generated by recursive parameter routines. 

We are now going to consider another requirement of our recursive routine A, 
which will lead us to the notion of isoparametricity of A. 

Let us turn back to the previous situation at the node p of the input circuit 
[3. Notations and assumptions will be the same as before. From Lemma [4] we 
deduce that the intermediate result of (3 associated with the node p, say G p , is 
a polynomial in X\, . . . ,X n whose coefficients form the entries of a geometrically 
robust, constructible map defined on M., say 9 p . Let T p be the image of this map and 
observe that T p is a constructible subset of a suitable afline space. The intermediate 
results of the circuit A((3) at the elements of the hypernode p of %a((3) constitute a 
polynomial vector which we denote by F p . 

We shall now make another requirement on the routine A at the node p of (3. 
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(B) There exists a geometrically robust, constructible map a p defined on T p such 
that a p o 6 P constitutes the coefficient vector of F p . 

In view of the comments made in |HKR11] . Section 3.3.1 we call the recursive 
routine A isoparametric (on input 0) if requirements (A) and (B) are satisfied at 
any node p of /3 which depends at least on one input. 

Suppose again that the recursive routine A is well behaved under restrictions. 
We call A well behaved under reductions (on input (3) if A(/3) satisfies the following 
requirement: 

Let p and p' be distinct nodes of (3 which compute the same intermediate 
results. Then the intermediate results at the hypernodes p and p' ofH.M$) 
are identical. Mutatis mutandis the same is true for the computation of 
the parameters of A at any node of ft. 

Assume that the routine A is recursive and well behaved under reductions. One 
verifies then easily that, taking into account the hypergraph structure Haw) °f ^-(P), 
any reduction procedure on /3 may canonically be extended to a reduction procedure 
ofA(p). 

It can be shown that under very light assumptions well behavedness under reduc- 
tions implies isoparametricity. Since from the point of view of software architecture 
well behavedness under reductions is a well motivated quality attribute of recursive 
routines we see that isoparametricity is computationally meaningful concept. For 
details we refer to |HKRllj . Section 3.3.2.3. 

5.2 Operations with recursive routines 

Let A and B be recursive routines as before and suppose that they are well behaved 
under restrictions and isoparametric or even well behaved under reductions. Assume 
that A({3) is an admissible input for B. We define the composed routine B o A in 
such a way that (B o A)((3) becomes the parameterized arithmetic circuit B(A(fl)). 
Since the routines A and B are well behaved under restrictions, we see easily that 
{B o A)((3) is a consistent, robust parameterized arithmetic circuit with parameter 
domain A4 . From Lemma [4] and Corollary [3] we deduce that B o A is a isoparametric 
recursive routine if A and B are isoparametric. In case that A and B are well 
behaved under reductions, one verifies immediately that B o A is also well behaved 
under reductions. Therefore, under these assumptions, we shall consider Bo A also 
SIS db routine of our computation model. 

The identity routine is trivially well behaved under restrictions and reductions 
and in particular isoparametric. 

Let A and B be two isoparametric recursive routines of our computation model. 
Assume that the robust parameterized arithmetic circuit is an admissible input 
for A and B and that there is given a one-to-one correspondence A which identifies 
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the output nodes of A((3) with the input nodes of B((3). Often, for a given input 
circuit (3, the correspondence A is clear by the context. If we limit ourselves to 
input circuits (3 where this occurs, we obtain from A and B a new routine, called 
their join, which transforms the input circuit (3 into the output circuit B((3) *\A{(3) 
(here we suppose that B(f3) *\ A((3) is consistent). Analyzing now B((3) *\ A((3), 
we see that the join of A with B is well behaved under restrictions in the most 
obvious sense. Since by assumption the routines A and B are recursive, the circuits 
A{(3) and B{(3) inherit from (3 a superstructure given by the hypergraphs T-La(/3) an d 
Hb{P)- Analyzing again this situation, we see that any reduction procedure on (3 
can be extended in a canonical way to the circuit B((3) *\ A((3). This means that 
the join of A with B is also well behaved under reductions if the same is true for 
A and B. More caution is at order with the notion of isoparametricity. In a strict 
sense, the join of two isoparametric recursive routines A and B is not necessarily 
isoparametric. However, condition (B) is still satisfied between the output nodes of 
(3 and B((3) *\ A{j3). A routine with this property is called output isoparametric. 

The union of the routines A and B assigns to the input circuit (3 the juxtaposition 
of A{(3) and B((3). Thus, on input (3, the final results of the union of A and B are the 
final results of A((3) and B{(3) (taken separately in case of ambiguity). The union 
of A and B behaves well under restrictions and reductions and is isoparametric if 
the same is true for A and B. 

Observe also that for a recursive routine A which behaves well under restrictions 
and reductions the following holds: let (3 be a robust parameterized arithmetic circuit 
that broadcasts to a circuit (3* and assume that (3 and (3* are admissible circuits for 
A. Then A{j3) broadcasts to A(f3*). 

5.3 Elementary routines 

From these considerations we conclude that routines, constructed as before by iter- 
ated applications of the operations isoparametric recursion, composition, join and 
union, are still, in a suitable sense, well behaved under restrictions and output 
isoparametric. If only recursive routines become involved that behave well under 
reductions, we may also allow broadcastings at the interface of two such operations. 

This remains true when we introduce, as we shall do now, in our computational 
model the following additional type of routine construction. 

Let (3 be the robust, parameterized circuit considered before, and let R(W; Y; X) 
be a generic computation belonging to our shape list. Let wp be a precomputed 
vector of geometrically robust constructible functions with domain of definition Ai 
and suppose that wp and W have the same vector length and that the entries of 
uifj are the final results of an output isoparametric parameter routine applied to the 
circuit (3. Moreover suppose that the final results of (3 form a vector of the same 
length as Y. 

Let K be the image of wp. Observe that K is a constructible subset of the afnne 
space which has the same dimension as the vector length of W. Denote by k the 
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vector of the restrictions to /C of the canonical projections of this affine space. We 
denote by R(k; Y; X) the ordinary arithmetic circuit over C obtained by substituting 
in the generic computation R(W;Y; X) the vector of parameter variables W by k. 
We shall now make the following requirement: 

(C) The ordinary arithmetic circuit R(k;Y; X) should be consistent and robust. 

Observe that requirement (C) is obsolete when R(W; Y; X) is a totally division- 
free ordinary arithmetic circuit. 

Suppose now that requirement (C) is satisfied. A new routine, say B, is obtained 
in the following way: on input ft the routine B joins with the generic computation 
R(W; Y\ X) at W and Y the previous computation of Wp and the circuit ft. 

From Lemma[4]and Corollary [3] we deduce that the resulting parameterized arith- 
metic circuit B(ft) has parameter domain hA and is robust if it is consistent. We 
shall therefore require that B(ft) is consistent (this condition is automatically satis- 
fied if R(k; Y\ X) is essentially division-free). One sees immediately that the routine 
B is well behaved under restrictions and reductions and is output isoparametric. 

From now on we shall always suppose that all our recursive routines are isopara- 
metric and that requirement (C) is satisfied when we apply this last type of routine 
construction. 

An elementary routine of our computation model is finally obtained by the 
iterated application of all these construction patterns, in particular the last one, 
isoparametric recursion, composition, join and union. As far as only recursion be- 
comes involved that is well behaved under reductions, we allow also broadcastings 
and reductions at the interface of two constructions. Of course, the identity routine 
belongs also to our model. The set of all these routines is therefore closed under 
these constructions and operations. 

We call an elementary routine essentially division-free if it admits as input only 
essentially division-free, robust parameterized arithmetic circuits and all specialized 
generic computations used to compose it are essentially division-free. The outputs 
of essentially division-free elementary routines are always essentially division-free 
robust circuits. The set of all essentially division-free elementary routines is also 
closed under the mentioned constructions and operations. 

We have seen that elementary routines are, in a suitable sense, well behaved 
under restrictions. In the following statement we formulate explicitly the property 
of an elementary routine to be output isoparametric. This will be fundamental in 
our subsequent complexity considerations. 

Proposition 7. Let A be an elementary routine of our computation model. Then A 
is output isoparametric. More explicitly, let ft be a robust, parameterized arithmetic 
circuit with parameter domain M.. Suppose that ft is an admissible input for A. Let 
9 be a geometrically robust, constructible map defined on A4 such that 9 represents 
the coefficient vector of the final results of ft and let T be the image of 9. Then T is a 
constructible subset of a suitable affine space and there exists a geometrically robust, 
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constructible map a defined on T such that the composition map a 06 represents the 
coefficient vector of the final results of A(/3). 

In case that A is a recursive routine, Proposition [7] expresses nothing but the 
requirement (B) applied to the output nodes of (3. 

5.4 End of the description of our computation model 

Elementary routines do not contain branchings. In order to capture the whole spec- 
trum of all known circuit based elimination methods in effective algebraic geometry, 
we are now going to introduce the concepts of algorithm and procedure admitting 
some limited branchings. This finishes the description of our computation model. 
We shall proceed rather informally. 

An algorithm will be a dynamic DAG of elementary routines which will be in- 
terpreted as pipes. At the end points of the pipes, decisions may be taken which 
depend on testing the validity of suitable universally quantified Boolean combina- 
tions of equalities between robust constructible functions defined on the parameter 
domain under consideration. The output of such an equality test is a bit vector 
which determines the next elementary routine (i.e., pipe) to be applied to the out- 
put circuit produced by the preceding elementary routine (pipe). This gives rise to a 
computation model which contains branchings. These branchings depend on a lim- 
ited type of decisions at the level of the underlying abstract data type, namely the 
mentioned equality tests. Because of this limitation of branchings, we shall call the 
algorithms of our model branching parsimonious (compare |GH01] and [CQH+03] ). 
A branching parsimonious algorithm A which accepts a robust parameterized arith- 
metic circuit with parameter domain Ai as input produces a new robust circuit 
A{0) with parameter domain Ai. In particular A(0) does not contain any branch- 
ings. 

Recall that our two main constructions of elementary routines depend on a pre- 
vious selection of generic computations from a given shape list. This selection may 
be handled by calculations with the indexing of the shape list. We shall think that 
these calculations become realized by deterministic Turing machines. At the begin- 
ning, for a given robust parametric input circuit with parameter domain Ai, a 
tuple of fixed (i.e., of (3 independent) length of natural numbers is determined. This 
tuple constitutes an initial configuration of a Turing machine computation which de- 
termines the generic computations of our shape list that intervene in the elementary 
routine under construction. The entries of this tuple of natural numbers are called 
invariants of the circuit 0. These invariants, whose values may also be Boolean (i.e., 
realized by the natural numbers or 1), depend mainly on algebraic or geometric 
properties of the final results of 0. However, they may also depend on structural 
properties of the labelled DAG 0. 

For example, the invariants of may express that has r parameters, n inputs 
and outputs, (over C) non-scalar size and depth at most L and /, that is totally 
division-free, that the final results of have degree at most d < 2 l and that for 



23 



any parameter instance their specializations form a reduced regular sequence in 
C[Xt, . . . , X n ], where X±, . . . , X n are the inputs of (3. 

Some of these invariants (e.g., the syntactical ones like number of parameters, 
inputs and outputs and non-scalar size and depth) may simply be read-off from the 
labelled DAG structure of f3. Others, like the truth value of the statement that the 
specializations of the final results of (3 at any parameter instance form a reduced 
regular sequence, have to be precomputed by an elimination algorithm from a previ- 
ously given software library in effective commutative algebra or algebraic geometry 
or their value has to be fixed in advance as a precondition for the elementary routine 
which becomes applied to f3. 

In the same vein we may equip any elementary routine A with a Turing com- 
putable function which from the values of the invariants of a given input circuit (3 
decides whether (3 is admissible for A, and, if this is the case, determines the generic 
computations of our shape list which intervene in the application of A to (3. 

We shall now go a step further letting depend the internal structure of the com- 
putation on the circuit (3. In the simplest case this means that we admit that the 
vector of invariants of (3, denoted by inv(/3), determines the architecture of a first 
elementary routine, say Ai nv (p), which admits (3 as input. Observe that the archi- 
tectures of the elementary routines of our computation model may be characterized 
by tuples of fixed length of natural numbers. We consider this characterization as 
an indexing of the elementary routines of our computation model. We may now use 
this indexing in order to combine dynamically elementary routines by composition, 
join and union. Let us restrict our attention to the case of composition. In this 
case the output circuit of one elementary routine is the input for the next routine. 
The elementary routines which compose this display become implemented as pipes 
which start with a robust input circuit and end with a robust output circuit. Given 
such a pipe and an input circuit 7 for the elementary routine B representing the 
pipe, we may apply suitable equality tests to the final results of -6(7) in order to 
determine a bit vector which we use to compute the index of the next elementary 
routine (seen as a new pipe) which will be applied to B(j) as input. 

A low level program of our extended computation model is now a text, namely 
the transition table of a deterministic Turing machine, which computes a function 
■if) realizing the following tasks. 

Let as before (3 be a robust parameterized arithmetic circuit. Then ifj returns 
first on input inv(/3) a Boolean value, zero or one, where one is interpreted as the 
informal statement u j3 is an admissible input". If this is the case, then tp returns 
on inv(/3) the index of an elementary routine, say A mv (p), which admits (3 as input. 
Then ip determines the equality tests which have to be realized with the final results 
of ^4inv( / 3)(/9)- Depending on the outcome of these equality tests ip determines an 
index value corresponding to a new elementary routine which admits Anv(/3)(/9) as 
input. Continuing in this way one obtains as end result an elementary routine 
A^\ which applied to (3, produces a final output circuit A^((3). The function ip 
represents all these index computations. We denote by ip(/3) the dynamic vector of 
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all data computed by ip on input (3. 

The algorithm represented by if) is the partial map between robust parametric 
arithmetic circuits that assigns to each admissible input (3 the circuit A>^(/3) as 
output. Observe that elementary routines are particular algorithms. This kind of 
algorithms constitute our computation model. We remark that any algorithm of this 
model is output isoparametric. If the pipes of an algorithm are all represented by 
essentially division-free elementary routines, we call the algorithm itself essentially 
division-free. 

One sees easily that the "Kronecker algorithm" |GLS01j (compare also |GHM + 98] , 
[GH H+97] and |GHMP97p for solving non-degenerate polynomial equation systems 
over the complex numbers may be programmed in our extended computation model. 
Observe that the Kronecker algorithm requires more than a single elementary rou- 
tine for its design. In order to understand this, recall that the Kronecker algorithm 
accepts as input an ordinary division-free arithmetic circuit which represents by 
its output nodes a reduced regular sequence of polynomials G\, . . . , G n belonging to 
C[X[, . . . , X n ]. In their turn, the polynomials G\, . . . , G n determine a degree pattern, 
say A := (Si, . . . , S n ), with Si := deg{Gi = 0, . . . , Gi — 0} for 1 < i < n. 

After putting the variables X\, . . . , X n in generic position with respect to G\, . . . , 
G n , the algorithm performs n recursive steps to eliminate them, one after the other. 
Finally the Kronecker algorithm produces an ordinary arithmetic circuit which com- 
putes the coefficients of n + 1 univariate polynomials P,V\, . . . ,V n over C. These 
polynomials constitute a "geometric solution" (see [GLS01] ) of the equation system 
G\ = 0, . . . , G n = because they represent the zero dimensional algebraic variety 
V := {Gi = 0, . . . , G n = 0} in the following "parameterized" form: 

V:={(V 1 (t),...,V n (t));teC,P(t) = 0}. 

Let (3 be any robust, parameterized arithmetic circuit with the same number of 
inputs and outputs, say X%, . . . , X n and G%(U, X\, . . . , X n ), . . . , G n {U, X%, . . . , X n ), 
respectively. Suppose that the parameter domain of (3, say M., is irreducible and 
that inv(/3) expresses that for each parameter instance u G Ai the polynomi- 
als G\{u, Xi, . . . , X n ), . . . , G n {u, X\, . . . , X n ) form a reduced regular sequence in 
C[Xi, . . . ,X n ] with fixed (i.e., from u G M. independent) degree pattern. Suppose, 
furthermore, that the degrees of the individual polynomials G\{u,X\, . . . , X n ), . . . , 
G n (u,Xi, . . . ,X n ) are also fixed and that the variables X±, . . . ,X n are in generic 
position with respect to the varieties {G\{u, X) = 0, . . . ,Gi{u,X) = 0}, 1 < i < n. 
Then, on input (3, the Kronecker algorithm runs a certain number (which depends on 
A) of elementary routines of our computation model which finally become combined 
by consistent iterative joins until the desired output is produced. 

We say that a given algorithm A of our extended model computes (only) param- 
eters if A satisfies the following condition: 

for any admissible input (3 the final results of A((3) are all parameters. 
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Suppose that A is such an algorithm and (3 is the robust parametric arithmetic 
circuit with parameter domain Ai which we have considered before. Observe that 
A(f3) contains the input variables X 1 , . . . , X n and that possibly new variables, which 
we call auxiliary, become introduced during the execution of the algorithm A on in- 
put (3. Since the algorithm A computes only parameters, the input and auxiliary 
variables become finally eliminated by the application of recursive parameter rou- 
tines and evaluations. We may therefore collect garbage in order to reduce A((3) to 
a final output circuit v4g na i(/3) whose intermediate results are only parameters. 

If we consider the algorithm A as a partial map which assigns to each admissible 
input circuit (3 its final output circuit v4.fi na i(/3), we call A a procedure. 

In this case, if ip is a low level program defining A, we call tp a low level procedure 
program. 

A particular feature of our extended computation model is the following: 
there exists a non-negative integer / (depending on the recursion depth of ^4) and 
non-decreasing real valued functions C/>0,...,C o >0 depending on one and the 
same dynamic integer vector, such that with the previous notations and Lp, Lj^p) 
denoting the non-scalar sizes of the circuits f3 and A(/3) the condition 

L AW) <C f W)H + --- + C ^(p)) 

is satisfied. 

In the case of the Kronecker algorithm (and most other elimination algorithms 
of effective algebraic geometry) we have / := 1, because the recursion depth of the 
basic routines which intervene is one. 

5.5 Procedures 

In the sequel we shall need a particular variant of the notion of a procedure 
which enables us to capture the following situation. 

Suppose we have to find a computational solution for a formally specified gen- 
eral algorithmic problem and that the formulation of the problem depends on cer- 
tain parameter variables, say Lq, . . . , U r , input variables, say Xi, . . . , X n and output 
variables, say Y 1 , . . . , Y s . Let such a problem formulation be given and suppose that 
its input is implemented by the robust parameterized arithmetic circuit (3 consid- 
ered before, interpreting the parameter variables Lq, . . . , U r as the basic parameters 

7Ti, . . . , 7T r . 

Then an algorithm A of our extended computation model which solves the given 
algorithmic problem should satisfy the architectural requirement we are going to 
describe now. 

The algorithm A should be the composition of two subalgorithms A^ and A^ 
of our computation model which satisfy on input (3 the following conditions: 

(i) The subalgorithm A^ computes only parameters, (3 is admissible for A^ and 
none of the indeterminates Y ± , . . . ,Y S is introduced in A^((3) as auxiliary 
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variable (all other auxiliary variables become eliminated during the execution 
of the subalgorithm A^ 1 ' on the input circuit (3). 

(ii) The circuit Afij] al ((3) is an admissible input for the subalgorithm A^ 2 \ the in- 
determinates K 1; . . . , Y s occur as auxiliary variables in A^ 2 \A^ al (f3)) and the 
final results of A^ 2 \A^l al (f3)) depend only on 7i"i, . . . , 7r r and Y l7 . . . , Y s . 

To the circuit A^ (A^^^)) we may, as in the case when we compute only pa- 
rameters, apply garbage collection. In this manner A^ (A^aiW)) becomes reduced 
to a final output circuit Aa na \(/3) with parameter domain M. which contains only 
the inputs Yi, . . . , Y s . 

Observe that the subalgorithm A^ is by Proposition [7] an output isoparamet- 
ric procedure of our extended computation model (the same is also true for the 
subalgorithm A^ 2 \ but this will not be relevant in the sequel). 

We consider the algorithm A, as well as the subalgorithms A^ and A^ 2 \ as 
procedures of our extended computation model. In case that the subprocedures 
A^ and A^ are essentially division-free, we call also the procedure A essentially 
division-free. 

The architectural requirement given by conditions (i) and (ii) may be interpreted 
as follows: 

the subprocedure A^ is a pipeline which transmits only parameters to the sub- 
procedure A^ 2 \ In particular, no (true) polynomial is transmitted from .A^ 1 ) to 
A^. 

Nevertheless, let us observe that on input (3 the procedure A establishes by 
means of the underlying low level program ip an additional link between (3 and the 
subprocedure A^ applied to the input A^((3). The elementary routines which 
constitute A^ on input A^ l \(3) become determined by index computations which 
realizes ip on inv(/3) and certain equality tests between the intermediate results of 
A {1) {(3). In this sense the subprocedure A^ transmits not only parameters to the 
subprocedure but also a limited amount of digital information which stems from the 
input circuit j3. 

The decomposition of the procedure A into two subprocedures A^ and A^ 
satisfying conditions (i) and (ii) represents an architectural restriction which is 
justified when it makes sense to require that on input the number of essential 
additions and multiplications contained in v4fi na i(/3) is bounded by a function which 
depends only on inv(/3). This is the case in the example of the next subsection where 
a substantial use of this restriction is made (see |HKRllj . Section 4.1, Observations). 

5.6 A hard elimination problem 

Let n e N and Si, ... , S„, T, Ui, . . . ,U n and Xi, . . . , X n be indeterminates. Let 
U := (U lt ...,U n ), S := (Si,...,S n ), X := (X 1 ,...,X n ) and := X\-X x - 
Si, ... , := X 2 n -X n - S n , ffW := £ 1<i<n 2^X t + T \\ l<l<n (\ + (17, - \)X % ). 
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Observe that the polynomials G\ , . . . , Gn form a reduced regular sequence in 
C[S, T, [/, X] and that they define a subvariety 14 of the affine space A™ x A 1 x A n x 
A n which is isomorphic to A™ x A 1 x A n and hence irreducible and of dimension 
2n + 1. Moreover, the morphism V n — > A n x A 1 x A n which associates to any point 
(s, t, u, x) e V n the point (s, t, u), is finite and generically unramified. Therefore the 
morphism ir n : V n — > A n x A 1 x A n x A 1 which associates to any (s, t, u, x) € V n 
the point (s,t,u,H^(t,u,x)) 6 A" x A 1 x A" x A 1 is finite and its image ir n (V n ) 
is a hypersurface of A n x A 1 x A™ x A 1 with irreducible minimal equation £ 
C[S, T,U,Y}. 

Thus, is an irreducible elimination polynomial of degree 2 n . Therefore any 
equation of C[S, T, U, Y] which defines 7i n (Vn) m A n x A 1 x A n x A 1 is up to a scalar 
factor a power of F^ n \ 

The equations G^f 1 = 0, . . . , G^ = and the polynomial represent a so 
called flat family of zero-dimensional elimination problems with associated elimi- 
nation polynomial (see [H KR1 1] . Section 4.1 for the notion of a flat family of 
zero-dimensional elimination problems). 

We consider now S±, . . . , S n , T, Ui, . . . , U n as basic parameters, X±, . . . , X n as 
input and Y as output variables. 

Let A be an essentially division-free procedure of our extended computation 
model satisfying the following condition: 

A accepts as input any robust parameterized arithmetic circuit which represents 
a flat family of zero-dimensional elimination problem with associated elimination 
polynomial F and A& n a\((3) has a single input Y and a single final result which 
defines the same hypersurface as F. 

With this notions and notations we have the following result. 

Theorem 8. There exist an ordinary division-free arithmetic circuit (3 n of size 
0(n) over C with inputs Si, ... , S n , T , Ui, . . . , U n , Xi, . . . , X n and final results 
G^, . . . , G^n\ H^ n \ The essentially division-free, robust parameterized arithmetic 
circuit 7„ := Afi na i((3 n ) depends on the basic parameters Si, ... , S n , T , U\, . . . ,U n 
and the input Y and its single final result is a power of F^ . The circuit 7 n performs 
at least fl(2 2 ) essential multiplications and at least Q(2 n ) multiplications with param- 
eters. As ordinary arithmetic circuit over C with inputs Si, ... , S n , T , U\, . . . ,U n 
and Y , the circuit 7 n has non-scalar size at least f2(2 n ). 

The proof of Theorem [8] is similar as that of Theorem [5] Moreover, Theorem [8] 
implies the asymptotic optimality of the Kronecker algorithm within our computa- 
tion model. For details we refer the reader to |HKRllj . Section 4 where also other 
examples of elimination problems are exhibited which are hard for algorithms of our 
computation model. 
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6 Approximative computations 



We are now going to apply the algorithmic model of Section [5] in a different context, 
namely that of approximative computations. 

Let (3 be a robust parameterized arithmetic circuit with parameter domain Ai, 
basic parameters 7Ti, . . . , n r and inputs X±, . . . , X n . Let U\, . . . , U r be parameter 
variables, U := (U%, . . . , U r ), n := (tti, . . . , ir r ), X := (Xi, . . . , X n ) and suppose that 
(3 is essentially division-free and has a single final result G. 

In this section we are going to introduce the notion of an approximative 0- 
computation and to discuss how in our computation model an approximative /3- 
computation can be transformed in an exact one. 

Let o be the vanishing ideal of Ai in C[U] and let us fix a polynomial P £ C[U] 
such that AAp is a Zariski open and dense subset of Ai . Let e be a new indeterminate. 

Definition 6 (Approximative parameter instance). An approximative parameter 
instance for (3 is a vector u(e) = («i(e), . . . , u r (e)) G C((e)) r which constitutes a 
meromorphic map germ at the origin such that any polynomial of a vanishes at u(e) 
and P(u(e)) ^ holds. 

Let u(e) be an approximative parameter instance for /3. Then there exists an 
open disc A around such that for any complex number c G A — {0} the germ 
u(e) is holomorphic at c and such that P(u(c)) ^ holds. This implies that any 
polynomial of a vanishes at u(c) and that in particular u(c) belongs to Ai. 

For technical reasons we need the following result. 

Lemma 9. Let notations and assumptions be as before. Let cf) : Ai — > A m be a 
geometrically robust constructible map and let u(e) be an approximative parameter 
instance for (3. Then there exists an open disc AofC around the origin and a germ 
ip of meromorphic functions at the origin such that u(e) and ip are holomorphic on 
A — {0} and such that any complex number c G A — {0} satisfies the conditions 
P(u{c)) ^ and ^(c) = <j>{u{c)). 

Proof. There exists an open disc A' of C around the origin such that u(e) is every- 
where defined on A' — {0} and such that any c G A' — {0} satisfies the condition 
P(u(c)) 7^ 0. Let Af be the Zariski closure of the image of A' — {0} under u(e). 
There exists a Zariski open and dense subset U of Af with U C Ai such that <fi is 
rational and everywhere defined on IA. Moreover there exists a non-zero polyno- 
mial Q G C[C/] such that Mq is contained in U and Zariski dense in Af. Therefore 
there exists a complex number cq G A' — {0} with Q(u(cq)) ^ 0. Hence the set 
K := {c G A' — {0};P(w(c)) = 0} is finite. Thus we may chose an open disc A 
around the origin with A C A' and A D K = such that the image of A under u(e) 
is contained in Mq. We conclude now that u(e) is everywhere defined on A — {0} 
and that every c G A — {0} satisfies the condition u(c) G U. For c G A — {0} let 
?/>(c) := (f){u{cj). Then ip : A — {0} — > C m is a well defined meromorphic function. 
Let c G A — {0} and let (ck)kem be a sequence of complex numbers c\. G A — {0} 
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which converges to c. Then (u(ck))ken is a sequence of points u(ck) G U which 
converges to u(c) and hence the sequence (<p(u(ck)))keN converges to <f>(u(c)) and is 
therefore bounded. This implies that ip is holomorphic at c. Thus ip is holomorphic 
on A — {0} and for any c G A — {0} we have ip(c) = <p(u(c)). Therefore we may 
interpret ip as a meromorphic map germ at the origin. □ 

Let u(e) be an approximative parameter instance. Then following Lemma [9] 
there exists an open disc A of C around the origin such that for any node p of (3 
with intermediate result G p (n,X) the expression G p (u(e),X) defines a polynomial 
in Xi, . . . ,X n whose coefficients are well determined meromorphic functions on A 
which are holomorphic on A — {0}. We denote by (3^ u ^ the labelled DAG of (3 
where we assign to each node p of (3 the polynomial G p (u(e),X). We call /?w e )) an 
approximative (3 -computation and denote by G^ u ^ the final result of (3^ u ^\ 

We say that the approximative /3-computation (3^ u ^ represents the polynomial 
H G C[X] if there exists a polynomial H^ u ^ G C[[e]][X] whose coefficient vector 
with respect to X constitutes a germ of functions which are holomorphic at the origin 
such that the final result of /3 (u(e)) can be written as G {u{e)) = H + e# (u(e)) . 

Let Wp be the set of coefficient vectors of the final results of the ordinary arith- 
metic circuits (3^ in C[X], where u G M. We consider Wp as a constructible subset 
of a suitable affine ambient space. Consequently the (Zariski or strong) closure Wp 
of Wp in its ambient space is well defined. In view of the next result the expression 
"approximative /3-computation" becomes selfexplanatory. 

Theorem 10. Let notations and assumptions be as before and let H G C[X]. Then 
the following three conditions are equivalent 

{%) there exists an approximative (3 -computation that represents H . 

(ii) there exists a sequence (uk)kem with Uk G M. such that the final results of the 
sequence ((3^ Uk ^)k^N of ordinary circuits converge to H in C[X]. 

(Hi) the coefficient vector of H belongs to Wp. 

Proof. The conditions (ii) and (Hi) are obviously equivalent because one is only a 
restatement of the other. It suffices therefore to show the implications (i) =^> (ii) 
and (in) =^ (i). We first prove (i) =^ (ii). 

Suppose that there exists an approximative parameter instance u(e) for (3 such 
that (3^ u ^ represents H G C[X] by means of a polynomial H^ u ^ G C[[e]][X] whose 
coefficient vector constitutes a holomorphic map germ at the origin. Thus H + eH u ^ 
is the final result of (3^ u ^\ We may choose a sequence (ek)keN of non-zero complex 
numbers such that for any k G N the germ u(e) is holomorphic at and satisfies 
the condition P(u(ek)) ^ and such that (ek)km converges to zero. Without loss 
of generality we may suppose that the coefficients of H^ u ^ are holomorphic at tk 
for any fcGN. 
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Let Uk ■= w(efc). Then Uk belongs to M. and H + ekH^ e '\ek,X) is the fi- 
nal result of the ordinary arithmetic circuit /3^ Uk \ Moreover the sequence (H + 
ekH u ^(ek, X)) fcgN converges to H. Therefore condition [ii) is satisfied. 

The other implication is more cumbersome. We adapt the argumentation of 
|Ald84j (see also |Lic90| . §A) to our computation model. Let us now prove (Hi) =>• 

(')• _ 

Without loss of generality we may suppose that M. is an irreducible afline variety. 
Let 9 be the geometrically robust constructible map which assigns to each u E M. 
the coefficient vector of the final result of the circuit (3^ u \ Observe that 9(Aip) 

is Zariski dense in W p. Thus B := W p — 9(A4p) is a proper Zariski closed subset 
of Wp. Let q be the dimension of the irreducible affine variety Wp and let h be the 
coefficient vector of the polynomial H. By assumption we have h G W p. In case of 
q = we conclude h G Wp and we are finished. Therefore we may suppose without 
loss of generality q > 0. 

By Noether's Normalization Lemma there exists a surjective finite morphism of 
irreducible affine varieties A : Wp — > A q . Since B is a proper Zariski closed subset 
of Wp we have X(B) ^ A q . We may therefore choose a point z G A q — X(B). Let 
L be a straight line of A r which passes through X(h) and z. Then X(h) belongs 
to L and X(B) R L is a finite set. Since the morphism A is finite, the irreducible 
components of A _1 (L) are all closed subcurves of Wp which become mapped onto 
L by A. In particular there exists an irreducible component C of A _1 (L) which 
contains h. Since X(B) R L is finite we have C ^ B. Therefore C R B is finite too. 
Suppose that 9(M. P ) R C is finite. Then we may conclude that C R B is infinite, a 
contradiction. Therefore 9(A4p) RC is infinite. Hence the Zariski closure of 9~ 1 (C) 
in A4 contains an irreducible component V such that 9(Vp) is Zariski dense in C. 
Let q* := dimV and u G V P . Observe q* > and that B* := 9~ 1 (9(u)) R V is a 
proper Zariski closed subset of V. Again by Noether's Normalization Lemma there 
exists a surjective finite morphism of irreducible affine varieties A* : V — > A 9 *. 

Since B* is a proper Zariski closed subset of V we have X*(B*) ^ A 9 *. Therefore 
we may choose again a point z* G A 9 * — X*(B*) and a straight line L* of A 9 * which 
passes through X*(u) and z*. Thus X*(B*) R L* is a finite set and X*(u) belongs to 
L*. The irreducible components of (A*) _1 (L*) are all closed subcurves of V which 
become mapped onto L* by A*. In particular there exists an irreducible component 
C* of (A*) _1 (L*) which contains u. Since X*(B*) RL* is finite we conclude C* ^ B*. 
Moreover from C* irreducible, u G C* and u G Vp we conclude that (C*)p is Zariski 
dense in C* . Hence C* B* implies that there exists a point u* G {C*)p — B* . 
Thus we have 9(u*) ^ 9(u). Therefore 9((C*)p) is Zariski dense in C. 

In this way we have found two irreducible closed subcurves C* of Ai and C of 
Wp with [C*)p nonempty such that 9 maps (C*)p into C and 6 , ((C*)p) is Zariski 
dense in C. Moreover we have h G C. 

The restriction of 9 to (C*)p is again a geometrically robust constructible map. 
Therefore 9 induces a finite field extension C(C) C C(C*). 

Let 7?. be the integral closure of C[C] in C(C*). Then m may be extended to a 
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maximal ideal 9JT of 1Z. Observe that 1Z is finitely generated over C. Thus 1Z is the 
coordinate ring of a non-singular curve and 9JT defines a point of this curve. Thus 
TZyji is a regular C-algebra of dimension one and therefore there exists an embedding 
of TZrgi into the power series ring C[[e]] which maps any generator of the maximal 
ideal of TZ^n onto a power series of order one. Moreover the elements of TZ<xn become 
mapped onto power series which constitute holomorphic function germs at the origin. 
Hence the coordinate functions of C*, given by the restrictions of 71*1, ...,7r r to 
C*, can be represented by Laurent series %(e), . . . , u r (e) of C((e)) which constitute 
meromorphic function germs at the origin. Let u(e) := (wi(e), ■ ■ • ,u r (e)). Then we 
have P(u(e)) 7^ and 9(u(ej) is by Lemmata well defined meromorphic map germ 
at the origin. Moreover 9(u(e)) belongs to C[C] and hence to IZm- This implies 
that the entries of the vector 6(u(e)) are power series of C[[e]] which constitute 
holomorphic function germs at the origin. Moreover h — 9(u(e)) belongs to the 
maximal ideal of C[C] m and hence to that of TZ<m- This means that e divides the 
entries of h — 9(u(e)) in C[[e]]. Observe now that 9{u(e)) is the coefficient vector 
of the final result G^ of (3^1 We conclude now that there exists a polynomial 
//"(«(«)) g C[[e]][X], whose coefficients constitute holomorphic function germs at the 
origin, such that H + eH^ = G^ u ^ holds. Since C* is contained in Ai we have 
finally A(u(e)) = for any polynomial A £ a. □ 



Theorem 10 is the technical main result of |Ald84j . where an analogous statement 
is proved in the particular case M. := A r (see also |Lic90j . §A). 



In the sequel we retake the notations of the proof of Theorem 10 



Let u(e) be an approximative parameter instance for /3 such that (3^^ rep- 
resents a polynomial H G C[X] whose coefficient vector h belongs to Wp. Then, 
using interpolation, one may first compute the coefficient vector 9(u(e)) of the final 
result G( u ( e )) of f3 u ( e > and then recompute G^ u ^ from the entries of 9(u(e)) and the 
indeterminates X\, . . . , X n . We may compute H in the same way from the entries 
of h = lim^o 9(u(e)) and X\, . . . , X n . This algorithm may be interpreted as two es- 
sentially division-free procedures ^4^^ and A^ 1 of our extended computation model 
where becomes applied to the input circuit (3^^ and returns the parameter 
vector #(«(€)) and A^ becomes applied to the input ^(^(e)) = A^(^ w ^ e '') and 
returns Thus A {1) computes on input (3^^ only parameters. Therefore 

A := (A^\ A^) may be interpreted as an essentially division-free procedure of our 
extended computation model which accepts (3^ 6 ^ as input and returns G^ u ^ as 
output (observe that Xi, . . . ,X n play simultaneously the role of input and output 
variables). Suppose now more generally that there is given an essentially division- 
free procedure A = (A^\ A^) of our extended computation model which accepts 
(3 as input and returns as output a robust parameterized arithmetic circuit A& na i((3) 
whose single final result is the polynomial G. 

Like in Section [5] we consider now the circuit 7 := An nSL i((3) and the geomet- 
rically robust constructible map v := A^\(3) whose domain of definition is A4. 
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Observe that 7 is an essentially division-free, robust arithmetic circuit with pa- 
rameter domain Ai. Let S be the image of the geometrically robust constructible 
map v. Then S is a constructible subset of a suitable affine space A p and there 
exists a geometrically robust constructible map ip : S — > A m and a vector of m- 
variate polynomials u* such that for u* := ip o v the geometrically robust con- 
structible map u*(b>*) = u* o v* is the coefficient vector 9 of the final result G of 
P. Observe that S* := ip(S) := v*(M) is a constructible subset of A m and that 
uj*(S*) = ui*(u*(M)) = 9(M) = Wp holds. Thus we may interpret u* as a surjec- 
tive polynomial map lu* : S* -»■ VFg. In the terminology of jCGH+03j . Section 3.2, 
the constructible set S* is a data structure, Wp an (abstract) object class and w* 
a holomorphic encoding of Wyj by <S*, where Wp represents the set of polynomials 
{G {u) ;u e M}. Moreover m is called the size of the data structure S. In an analo- 
gous way we may also consider S as a data structure and u := u* o ijj an encoding 
of (observe that Wp = 9(M) = u* o ip o z/(A^) = u{S) holds). Again in the 
terminology of [C GH + 03~] . Section 3.2, a; is a continuous encoding of the object class 
Wp by the data structure S whose size is p. 

Lemma 11. Let (w fc ) fceN; u k G M, (s k ) km , s k G S and (s* k ) k( z N , s* k G S* three 
sequences which satisfy the following condition: 

the sequences (8(uk))keti> ( w ( s fc))fceN and (cu* (s* k )) k( z^ converge each to a 
point of Wp . 

Then the sequences (^(u k )) k ^, (s k ) k( z^ and (i^*(u k )) k( z^, (s* k ) k< =^ converge each to a 
point of S or S* , respectively. 

Proof. Observe that the procedure is output isoparametric. Therefore there 
exists a geometrically robust constructible map o with domain of definition Wp 
such that a o 9 = v holds. Since a and 9 are continuous with respect to the strong 
topologies of their source and image spaces and (9(uk))k& converges to a point g 
of Wp, the sequence (is(u k )) k £N converges to cr(g). Hence there exists a parameter 
instance mgM with g = 9(u). Thus we have o~(g) — (a o 9)(u) = u(u) and u(u) 
belongs to S. This implies that the sequence (v(u k )) ke ^ converges to a point of S. 

Observe that for each k G N there exists a parameter instance u k G M. with 
s k = v{u k ). Therefore we have 

u(s k ) = u) o u(u k ) = (u* o if; o u){u k ) = (u* o v*)(u k ) = 9(u k ). 

By assumption the sequence (cu(s k )) k£ ?>}, and hence also the sequence (9(u k )) ke ^, 
converge to a point of Wp. As we have shown, this implies that the sequence 
(z/(wfc))fc e N, i-e. the sequence (sfc)fc £ N converges to a point of S, namely to cr(g). 

The statements about the sequences (v*(u k )) ke ^ and (s£) fcgN are shown in the 
same way. □ 
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Lemma 11 implies that u and u* are robust encodings in the terminology of 
|CGH+03j . Section 3.2. 

Let now be given an approximative parameter instance u(e) of (3 such that /3W e )) 
represents a polynomial H whose coefficient vector belongs to Wp. 

We wish to evaluate the polynomial H. There exists a parameter instance u E Ai 
with H = G( u \ but we are in general not able to find it. However we may be able 
to find a sequence (ek)km of non zero complex numbers such that for any G N 
the germ u(e) is holomorphic at e k and satisfies the condition P(u(e^)) ^ 0. This 
implies that u k := u(e k ) belongs to M.. Since the approximative /3-computation 
/jOO)) represents the polynomial H, the sequence (G^^jteN converges to H. From 
Lemma 11 we deduce now that the sequences (v(uk))kem and (^*(wfc))fceN converge 
to points s and s* of S and S*, respectively. Therefore co(s) and u*(s*) are equal 
and form the coefficient vector of H. In case that the essential parameters of the 
circuit 7 = ^4fi na i(/3) constitute the entries of the geometrically robust constructible 
map u*, we may reinterpret 7 as a robust parameterized arithmetic circuit with 
parameter domain S*. In this interpretation j^ 8 ' becomes an ordinary division-free 
arithmetic circuit in C[X] whose single final result is H. 

If the coefficient vector of H belongs to Wp, but not to Wp we need the additional 
requirement that the sequences (v(u(ek)))km or ( Z7 *(^( e fc)))fceN are convergent (or at 
least bounded) in order to make the same kind of conclusions. However this depends 
on the technical condition that the geometrically robust constructible map a of the 



proof of Lemma 11 may continuously be extended to the closure Wp of Wp or on 
the choice of the sequence (e^ken- 

In order to illustrate these technical considerations we are now going to discuss 
the following example. 

Let L, n be natural numbers, r := (L + n + l) 2 , Kl jU '■= 4(L + n + l) 2 + 2 and 
let Ui, . . . ,U r be parameter and Xi, . . . , X n input variables. Let Aii,n '■= A r and 
7Ti, . . . , 7T r the canonical projections of A4l,u onto A 1 . 

Following [BC S97] , Chapter 9, Exercise 9.18 there exists a totally division-free 
generic computation /3^,n with a single final result Gl, h such that any polynomial 
H G C[Xi, . . . ,X n ] is evaluable by at most L essential multiplications if and only if 
there exists a parameter instance u G M.L,n such that H = G^/ n holds. As a totally 
division-free generic computation we may interpret (3l )71 as a robust parameterized 
arithmetic circuit with parameter domain M.L,n, basic parameters 71*1,..., 7iy and 
inputs Xi, . . . , X n . 

Let us fix points 71, ... , 7^ n of Z n such that (71, . . . , jk l n ) is an identification 
sequence for the set Wi, n '■= G M.L,n}- This means that for any Hi,H 2 G 

W L}n the implication #1(71) = H 2 (71 ),•••, #1(7^, J = H 2 {^K L>n ) =>- #1 = #2 
holds. 

There exist identification sequences (71, ... , 7^ n ) of bit length at most 4(L + 1) 
(see jCGH+03j and jGHMSllj for details). 

As before let A = (A^ l \A^) be an essentially division-free procedure of our 
computation model which on input 0L n returns a robust parameterized arithmetic 
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circuit Asnai(pL,n) whose single final result is Gl,u- Let Vl, u '■= A^ l \^L,n) an d let 
Sl,u be the image of the geometrically robust constructible map vl,u- We think 
Sl,u as a constructible subset of an affine space A PL ' n . Furthermore let be given a 
geometrically robust constructible map ij)L,n '■ <Sl,u A mi - n and a vector of m L>n - 
variate polynomials u* L n such that for v* L n := ipL,n° ^L,n the vector of coefficients of 
Gl,u with respect to the variables X\, . . . , X n can be written as u* L n o v* L n . 

Assume now pi, n = Kl,u and that z/£ )Tl : Aii,n — > <Sl,u is the polynomial (and 
hence geometrically robust constructible) map which for u G M.Ln is defined by 

^»(«):=(GW(7i),...,G (u) (7m £ ,J)- 

Let A^ 1 be a procedure of our computation model which accepts v L ^ n = A(^L >n ) 

as input and returns on u^ n as output the polynomial Gi, n (such procedures always 

exist). Then the size pi, n = Ki,n of the continuous encoding uj* Ln o n : Sl, u — > 

Wp L:n of W htn is K L<n = 4(L + n + l) 2 + 2 whereas for S* L>n := j L , n {S) = v* L JM) 

the map lo* l n : S* L n — >■ W 7 "^ n represents a holomorphic encoding of n of size 

^L,n = 2 n<yLn ' > (see [GHMSllJ, Theorem 23). The proof of this lower bound may be 

adapted to show the estimate = 2 n ( Ln * > for any procedure A = (A^^^) of 

our extended computation model which recalculates the polynomial Gl,u from the 

circuit pL,n- 

Thus, there are natural classes of polynomials which have continuous encodings 
of "small size" whereas their holomorphic encodings may become necessarily "large" . 



References 

[Ald84] A. Alder. Grenzrang und Grenzkomplexitat aus algebraischer und topol- 
ogischer Sicht. PhD thesis, Universitat Zurich, Philosophische Fakultat 
II, 1984. 

[BCS97] P. Biirgisser, M. Clausen, M. A. Shokrollahi. Algebraic Complexity The- 
ory. Grundlehren der mathematischen Wissenschaften, 315. Springer 
Verlag, 1997. 

[CGH89] L. Caniglia, A. Galligo, J. Heintz. Some new effectivity bounds in com- 
putational geometry. Applied Algebra, Algebraic Algorithms and Error 
Correcting Codes. Proc. of the 6th Intern. Conference AAECC, Best 
Paper Award AAECC-6. Springer LNCS, 357 131-151, 1989. 

[CGH+03] D. Castro, M. Giusti, J. Heintz, G. Matera, L. M. Pardo. The hardness 
of polynomial equation solving. Foundations of Computational Mathe- 
matics, 3 (4) 347-420, 2003. 

[DFGS91] A. Dickenstein, N. Fitchas, M. Giusti, C. Sessa. The membership prob- 
lem of unmixed polynomial ideals is solvable in single exponential time. 
Discrete Applied Mathematics, 33 73-94, 1991. 



35 



[FIK86] T. Freeman, G. Imirzian, E. Kaltofen. A system for manipulating poly- 
nomials given by straight-line programs. In Proceedings of the fifth ACM 
Symposium on Symbolic and Algebraic Computation, SYMSAC '86, 169- 
175. ACM, 1986. 

[GH01] M. Giusti, J. Heintz. Kronecker's smart, little black boxes. In Founda- 
tions of Computational Mathematics, R. A. DeVore, A. Iserles, E. Siili 
eds., 284 of London Mathematical Society Lecture Note Series, 69-104. 
Cambridge University Press, Cambridge, 2001. 

[GHH+97] M. Giusti, K. Hagele, J. Heintz, J. L. Montana, J. E. Morais, L. M. 

Pardo. Lower bounds for diophantine approximation. Journal of Pure 
and Applied Algebra, 117 277-317, 1997. 

[GHM+98] M. Giusti, J. Heintz, J.E. Morais, J. Morgenstern, L.M. Pardo. Straight- 
line programs in geometric elimination theory. Journal of Pure and 
Applied Algebra, 124 101-146, 1998. 

[GHMP95] M. Giusti, J. Heintz, J. E. Morais, L. M. Pardo. When polynomial 
equation systems can be "solved" fast? In Applied Algebra, Algebraic 
Algorithms and Error Correcting Codes. Proc. of the 11th Intern. Con- 
ference AAECC. LNCS, 948, 205-231. Springer, 1995. 

[GHMP97] M. Giusti, J. Heintz, J. E. Morais, L. M. Pardo. Le role des structures 
de donnees dans les problemes d'elimination. Comptes Rendus Acad. 
Sci., Serie 1 (325) 1223-1228, 1997. 

[GHMS11] N. Gimenez, J. Heintz, G. Matera, P. Solerno. Lower complexity bounds 
for interpolation algorithms. Journal of Complexity, 27 151-187, 2011. 

[GLS01] M. Giusti, G. Lecerf, B. Salvy. A Grobner Free Alternative for Polyno- 
mial System Solving. Journal of Complexity, 17 154-211, 2001. 

[Hei89] J. Heintz. On the computational complexity of polynomials and bilinear 
mappings. A survey. Proceedings 5th International Symposium on Ap- 
plied Algebra, Algebraic Algorithms and Error Correcting Codes Springer 
LNCS, 356 269-300, 1989. 

[HKR11] J. Heintz, B. Kuijpers, A. Rojas Paredes. Software engineering and 
complexity in effective algebraic geometry. Submited to J. of Com- 



plexity, |http : //arxiv . org/PS_cache/arxiv/pdf / 1 110/1110. 3030v2 . 
pfj 2011. 

[HM93] J. Heintz, J. Morgenstern. On the intrinsic complexity of elimination 
theory. Journal of Complexity, 9 471-498, 1993. 



36 



[HMW01] J. Heintz, G. Matera, A. Weissbein. On the time-space complexity of 
geometric elimination procedures. Applicable Algebra in Engineering. 
Communication and Computing (AAECC Journal), 11 239-296, 2001. 

[HS81] J. Heintz, M. Sieveking. Absolute primality of polynomials is decidable 
in random polynomial time in the number of variables. In S. Even, 
O. Kariv, editors, International Colloquium on Automata, Languages 
and Programming ICALP, 115 of Lecture Notes in Computer Science, 
16-28. Springer, 1981. 

[Ive73] B. Iversen. Generic Local Structure of the Morphisms in Commutative 
Algebra. Springer- Verlag, Berlin, 1973. 

[Kal88] E. Kaltofen. Greatest common divisors of polynomials given by straight- 
line programs. Journal of the ACM, 35 (1) 231-264, 1988. 

[Kro82] L. Kronecker. Grundziige Einer arithmetischen Theorie der algebrais- 
chen Grossen (Fundamentals of an arithmetic theory of algebraic quan- 
tities). J. Reine Angew. Math., 92 1-122, 1882. 

[Kun85] E. Kunz. Introduction to commutative algebra and algebraic geometry. 
Birkhauser, Boston, 1985. 

[Lan93] S. Lang. Algebra. Addison- Wesley, Massachusetts, 1993. 

[Lec] G. Lecerf. Kronecker: a Magma package for polynomial system 

solving, http : //lecerf . perso . math . cnrs . f r/ sof tware/kronecker/ 
inde x. htmll 

T. M. Lickteig. On semialgebraic decision complexity. Habilitationss- 
chrift, Universitat Tubingen TR-90-052, Int. Comp. Sc. Inst., Berkeley, 
1990. 

F. S. Macualay. Algebraic theory of modular systems. Tracts Math. 
Cambridge University Press, 1916. 

T. Mora. SPES I: The Kronecker- Duval Philosophy. Cambridge Uni- 
versity Press, 2003. 

T. Mora. SPES II: Macaulay's Paradigm and Groebner Technology. 
Cambridge University Press, 2005. 

D. Mumford. The red book of varieties and schemes, 1358. Springer, 
Berlin Heidelberg, New York, 1. edition, 1988. 

M. S. Paterson, L. J. Stockmeyer. On the number of nonscalar multipli- 
cations necessary to evaluate polynomials. SIAM Journal on Computing, 
2 60-66, 1973. 



[Lic90] 

[Macl6] 

[Mor03] 

[Mor05] 

[Mum88] 

[PS73] 



37 



[Sha94] I. R. Shafarevich. Basic algebraic geometry: Varieties in projective space. 
Springer, Berlin Heidelberg, New York, 1994. 

[SS95] M. Shub, S. Smale. On the intractability of Hilbert's Nullstellensatz and 
an algebraic version of "NP^P?". Duke Math. J., 81 47-54, 1995. 

[vdW50] B. L. van der Waerden. Modern Algebra II. Ungar, New York, 1950. 

[ZS60] O. Zariski, P. Samuel. Commutative algebra II, 39. Springer, New York, 
1960. 



38 



